Analysis of the Interface in a Nonequilibrium Two- Temperature Ising Model 
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The influence of nonequilibrium bulk conditions on the properties of the interfaces exhibited 
by a kinetic Ising-like model system with nonequilibrium steady states is studied. The system is 
maintained out of equilibrium by perturbing the familiar spin-flip dynamics at temperature T with 
completely-random flips; one may interpret these as ideally simulating some (dynamic) impuri- 
ties. We find evidence that, in the present case, the nonequilibrium mechanism adds to the basic 
thermal one resulting on a renormalization of microscopic parameters such as the probability of 
interfacial broken bonds. On this assumption, we develop theory for the nonequilibrium "surface 
tension", which happens to show a non-monotonous behavior with a maximum at some finite T. 
The phase diagram, as derived from this effective interfacial free energy, exhibits reentrant behavior. 
In addition, interface fluctuations differ qualitatively from the equilibrium case, e.g., the interface 
remains rough at zero-T, in full agreement with Monte Carlo simulations. We discuss on some 
consequences of these facts for nucleation theory, and make some explicit predictions concerning the 
nonequilibrium droplet structure. 
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I. INTRODUCTION 

Interfaces that separate different homogeneous media 
are familiar from many natural phenomena such as phase 
segregation, wetting processes, fluid dynamics, crystal 
growth, and molecular beam epitaxy, for instance. In 
practice, the interface may determine the system mor- 
phology and its critical properties, or the details of time 
evolution, which has motivated many specific studies dur- 
ing the last two decades. 0,1111 Ell lH0 The phenom- 
ena most deeply studied so far concern interfaces that 
separate equilibrium phases. However, the actual systems 
of interest are seldom at equilibrium and one often needs 
to be concerned with nonequilibrium interfaces, i.e., in- 
terfaces that separate nonequilibrium phases . [jl Hoi HH 
This paper aims towards a better understanding of how 
nonequilibrium conditions influence the properties of an 
interface. The results here will be applied in a forth- 
coming paper to analyze the exit from a (nonequilib- 
rium) metastable state, which is an interface-controlled 
process, [l^ 

Mathematical complexity often compels one to deal 
with the simplest model situations. In this paper, we 
study interfaces in a two-dimensional kinetic Ising model. 
The nonequilibrium condition is obtained by perturb- 
ing the underlying stochastic dynamics in such a way 
that, in general, a nonequilibrium steady state is reached 
asymptotically (instead of the more familiar thermody- 
namic equilibrium state). We develop for this case a 
simple approximation which predicts both microscopic 



and macroscopic properties for the (nonequilibrium) in- 
terface, namely, the profile or single-step height probabil- 
ity distribution and a (nonequilibrium) "surface tension" , 
<7 nc . The latter turns out to qualitatively differ from the 
equilibrium surface tension, o~ e . In particular, o~ ne behaves 
non-monotonously with decreasing temperature, and ex- 
hibits a maximum at some finite temperature, unlike o~ e 
that monotonously grows as one cools the system. In ad- 
dition, the phase diagram of the model, as derived from 
this effective surface tension, exhibits reentrant behavior. 
We also predict that, due to the nonequilibrium pertur- 
bation, the interface in this model remains rough at zero 
temperature. In order to understand this behavior, we 
analyze the shape of a droplet of the minority phase, and 
conclude how the nonequilibrium condition substantially 
influences the low-temperature droplet morphology. Our 
predictions are compared with the results from computer 
simulations, which seem to firmly support them in gen- 
eral. 

The paper is organized as follows. The model is defined 
in § II, and § III describes our approximation. The main 
results are in § IV which also contains some details of the 
computer simulations and a comparison of the numerical 
results with theory. § V is devoted to conclusions. 



II. THE MODEL 



Let the square lattice A(L x ,L y ) E Z 2 of size N = 
L x x L y with a binary spin variable, Si = ±1, at each node 
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i E [lj-^V]- We remind that the two states of Sj may be 
interpreted as corresponding to the presence or absence, 
respectively, of a particle at i; this happens to provide a 
more intuitive picture concerning the phenomena of inter- 
est here. There is interaction between nearest-neighbor 
spins given by the Ising Hamiltonian, 



H 



(i,3) 



(1) 



and stochastic dynamics by single-spin flips. The latter 
occur with transition rate (per unit time) given by 



Uj(Si 



p 



il-p)*(pATU). 



(2) 



Here, (3 — l/T — we take both the coupling constant and 
the Boltzmann constant equal to unity — , and AHi = 
4sj(ni — 2), where rii G [0,4] is the number of up 
nearest-neighbor spins surrounding s<, i.e., AHi mea- 
sures the energy cost of flipping at i. The undetermined 
function in (0) is either * (T) = e~ r (l + e~ r ) _1 or 
^ (r) = min [l,e~ r ] ; as indicated below, we are mostly 
concerned here with the first choice, except when the sec- 
ond one allows for an explicit or a simpler description. In 
any case, both choices lead to the results in this paper. 

One may interpret that the parameter p in J5J balances 
the competition between two thermal baths: one is at 
temperature T, while the other induces completely ran- 
dom transitions as if it was at infinite temperature. For 
p = 0, u(si — > —Si) satisfies detailed balance, and the 
system goes asymptotically to the equilibrium state for 
temperature T and energy H. The system exhibits in this 
case the familiar, Onsager critical point at T = Tcip = 
0) = To = 2/ln(l + V2) . Otherwise, < p < 1, the 
competition in @ impedes canonical equilibrium and, in 
general, a nonequilibrium steady state sets in asymptot- 
ically with time. A critical point is still observed in this 
case, but at T = Tc{p) < To, as far asp < p c . For p > p c 
the system remains in the disordered phase at any T. As 
shown below, our theoretical approach in this paper pre- 
dicts p c — (v2 — l) 2 ~ 0.1716, in perfect agreement with 
previous Monte Carlo estimations. I n addition, when 
subject to an external magnetic field, this model exhibits 
metastable states whose strength decreases with increas- 
ing field, eventually becoming unstable. The spinodal 
field characterizing the limit of metastability undergoes 
an interesting reentrant phenomenon as a consequence 
of the non- linear interplay between T and p. 01 We will 
show below that a similar reentrant behavior is observed 
for the phase diagram Tcip)- 

In the ordered phase below Tc{p), this system ex- 
hibits an interface for appropriate boundary conditions. 
Consider, for instance, periodic boundary conditions 
along the i-direction, and open boundaries along the 
y-dircction. Let then the spins in the bottom (top) 
row to be frozen in the up(down)-state, while the rest 
of spins are allowed to change stochastically according 
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FIG. 1: Example of interface for L x — 12 with steps S = 
(1, -2, -1, 2, 2, -2, 1, -2, 1, 1, -1, 0). The numbers shown in 
the squares indicate the class of the corresponding spin as 
defined in the main text. Notice that interfacial spins can 
only belong to classes 2, 3 and 4. 



to (J5J). Under these conditions, an interface eventu- 
ally develops along the ^-direction that separates up- 
from down-spin rich regions located at the bottom and 
top of the system, respectively. For p = (the equi- 
librium case), the macroscopic properties of this inter- 
face are well-known. In particular, its scaling behavior 
is characteristic of the Kardar-Parisi-Zhang universality 
class, and one knows the total free energy per unit 
length or surface tension. 0] The question is how these 
properties are affected by the nonequilibrium perturba- 
tion parametrized by p. 



III. SOLID-ON-SOLID APPROXIMATION 

The interface is first analyzed here by adapting to our 
case the solid-on-solid (SOS) picture introduced by Bur- 
ton, Cabrera and Frank. This assumes that the inter- 
face can be described by a single-valued discrete func- 
tion completely defined by the set of interface steps, 
{S x ,x S as illustrated in Fig. ^ No over- 

hangs are allowed in this approximation. Furthermore, 
the heights of the individual steps are assumed to be inde- 
pendent. The probability of a step of height S is assumed 
to be given by 



P(6) = 



1 



z[T,P, 7 (</>)] 



X{T,p)^^ Ws . 



(3) 



Here X(T,p) is the statistical weight associated to a bro- 
ken bond in the interface, and 7(</>) is a Lagrange mul- 
tiplier intended to keep the average step value at a x- 
independent value, (S x ) = tan0, where <fi is the aver- 
age angle between the interface and the x axis. 0, 0] 
The function z[T, p, 7] may be obtained from the nor- 
malization of and the "partition function" for the 



3 



SOS interface follows as Z(T,p,j) = [Xz(T,p, 7 )] L *. 
This allows one to define a "thermodynamic" potential 
ip(T,p,j) = limi x _ (00 — (fiLx)' 1 In Z. Our interest is 
then on the "free energy" a'(T,p,4>) — defined as con- 
jugated to the potential ip(T,p, 7) via a Legendre trans- 
form, which involves the variables tan</> and T7 — and, 
in particular, on its projection a along the x axis. The 
nature of the above approximation is discussed below; 
we now remark that our use of equilibrium words here is 
only for simplicity and comfort. The general result is 



surface tension Q may be now explicitly written as 



a((j)) = T J cos 



01 {7(0) 



tan 6— In 



X(l -A 2 



1 + X 2 - 2X cosh 7 (</>) J ' 
(4) 

0, one has the weight 



In the equilibrium limit p 

X(T,p = 0) = e" 2/3 , and a c = a(T,p = 0, < 
the SOS surface tension associated to the equilibrium 
Ising interface. [IH UtI Il8j| We shall assume in this paper 
that all the main effects of the nonequilibrium pertur- 
bation (p) on the interface can be taken into account 
after a proper generalization of the microscopic param- 
eter X(T,p). The function er no = o~(T,p > 0,0) which 
results after using this generalization in Q is therefore 
assumed to be the nonequilibrium "surface tension" in 
SOS approximation, and this is expected to capture the 
macroscopic properties of the nonequilibrium interface. 
The resulting cr ne is to be interpreted as an effective 
free-energy per unit length in the present case that lacks 
of a proper bulk free-energy function. This definition 
of nonequilibrium "surface tension" is based on the as- 
sumption that the normalization Z{T, p, 7) of the prob- 
ability measure associated to interface configurations in 
SOS picture is some sort of nonequilibrium analog of the 
partition function. Similar hypothesis have been shown 
to yield excellent results when applied to other nonequi- 
librium models. ,19| For instance, in the one-dimensional 
asymmetric simple exclusion process (ASEP) with open 
boundaries, the distribution of (complex) zeros of the 
steady-state normalization factor has been shown to 
obey the Lee- Yang picture of phase transitions. 

In the simplest scenario — which is consistent with the 
equilibrium limit — , the weight X(T,p) will not depend 
on 7(0). An explicit relation between 7(0) and tan0 can 
then be obtained. In particular, using (5 X ) = tan0, one 
finds that 



a±7(0) _ (l + A 2 )tan0±S(0) 



z(tf>) = 



2X(tan0± 1) 
(1 - X 2 )(l - tan 2 
1 + A 2 -S(0) 



(5) 
(0) 



where S(</>) = [(l-X 2 ) 2 tan 2 0+4A 2 ] 1 / 2 , and we dropped 
an obvious dependence on T and p. The nonequilibrium 



CT (</>) = T|cos0N (tan 0) In 



- In 



(1 + X 2 )tan0 + S*(0) 



2X(tan0 + 1) 



X(l- X 2 )(l-tan 2 
1 + X 2 - S(<j>) 



(7) 



Equations J2J and Q are two important properties of 
the interface at the microscopic and macroscopic levels 
of description, respectively. Interesting enough, re- 
duces in equilibrium (p = 0) to the known exact result 
for (f> = 0. It also yields a very good approximation for 
any angle |0| < 7r/4;0,|2(J for \(/)\ > n/4, it is convenient 
to turn to the y (instead of x) axis as the reference frame. 
o~(T,p = 0, 4> = 0) is a monotonously decreasing function 
of T, which converges toward 2 as T — > 0, and vanishes 
at the exact Onsager critical temperature. In addition, 
the angular dependence of cr(T,p — 0,0) may be used 
to determine the equilibrium crystal shape via minimiza- 
tion of the total surface tension for a fixed volume in a 
homogeneous droplet (Wulff construction). 



Statistical Weight of a Broken Bond 

At equilibrium, the weight equals the Boltzmann fac- 
tor, X(T,p = 0) — e~ 2/3 (2 is the energy cost of a broken 
bond). More generally, < p < 1, the weight of a broken 
bond depends on the local order surrounding the spin at 
the end of the bond. One may say that X(T,p) depends 
on the class of this spin. The spin Si is said to be of class 
ry(si,n;) = 3 — Si(ni — 2), where n.; is the number of up 
nearest-neighbor spins of Si. That is, our model may ex- 
hibit up to five different classes of spins, r\ = 1, . . . , 5.|2l) 
All the spins in a given class are characterized by the 
same value of AH V = 4si(n, — 2) and, consequently, by 
the same rate ©• 

The function X v (T, p) for class rj now follows straight- 
forwardly. Consider J2| with $ (r) = min [l,e _r ] which 
allows for a simpler and more explicit discussion. 22] If, 
for instance, a spin in class rj = 1 is flipped, four new 
broken bonds appear. Since ATCi = 8, one immediately 
has that Xi(T,p) = \p + (1 - p)e" 8 ' 3 ] 1 / 4 for the first 
class. Equivalents, X 2 {T,p) = \p + (1 - p)e" 413 ] 1 ^ 2 for 
the second class, and the rest are characterized by the 
same weight as in equilibrium, X^iT^p) = X^T^p) — 
X 5 (T,p) = X{T,p = 0). It follows that Ai < X 2 < 
^3A,5 f° r an y < p < 1, and one may also see that 
X v (T 7 p — * 0) — > e _2/3 , independent of 77, as expected. 

As illustrated in Fig. interfacial spins may only be- 
long to classes 2, 3 and 4; the class 1 corresponds to 
spins in the bulk of (either up or down) homogeneous 
regions, and the class 5 corresponds to typical isolated 
fluctuations in the bulk. Consequently, for the case in 
consideration,!^ only two weights, X 2 and A3 = A4 
are relevant. An even simpler description ensues assum- 



4 



Step variables 


Configuration 


P(S,e) x Q 


8 > 0, e > 


5 


J 

1 
1 


A S +tx S + e-2 x 3 


8 > 0, e = 


it 




(5 > 0, e < 


5 


E 


<5 + e ji^" A — a jj^" 2 a + 1 


5 = 0, e > 






5 = 0, e = 


5 e 




(5 = 0, e < 


8 


A e xi e|_1 X3 2 


(5 < 0, e > 




^<5 + ej£-A-aj^2a + l 


5<0, e = 




A Ij X? 


5 < 0, e < 


\ 


A «+e x |5|+|e|-2 x | 



TABLE I: The nine different typical configurations of an in- 
terfacial spin column in our approximation. These configu- 
rations are defined by the signs of the left, 5, and right, e, 
steps. The last column shows the probability P(8, e) of each 
configuration. Here, a = min(|5|,|e|) and A = max(|5|, |e|); 
see the main text for other definitions.. 



ing, which amounts a reasonable mean-field approxima- 
tion, that interfacial broken bonds have an unique statis- 
tical weight equal to the weighted average of X 2 and X 3 , 
namely, 

X(T,p)=IL 2 (T,p)X 2 (T,p) 

+ [Tl 3 (T,p) + Tl4T,p)]X 3 (T,p). (8) 

Here, H n (T, p) is the probability of an interfacial broken 
bond associated to a spin of class r\. Alternatively, given 
that any bond can be arbitrarily associated to any of the 
two spins at the ends, we may interpret n,,(T, p) as the 
probability of an interfacial up spin of class 77. Therefore, 
U 2 (T,p) + U 3 (T,p)+Ih(T,p) = l. 



Population of Interfacial Spin Classes 

Next, we estimate the population densities H v (T,p), 
which requires a detailed counting beyond the SOS ap- 
proximation. Let P(5, e) the joint probability that a step 
variable equals S and the step variable at its right is e. 



Step variables 


m(5, e) 


n 3 (5, e) 


114(5, e) 


N(S,e) 


8 > 0, e > 


5-1 


1 





8 


8 > 0, e = 


5-1 


1 





8 


8 > 0, e < 


A — a 


a - 1 


1 


A 


8 = 0, e > 


1 








1 


8 = 0, e = 


1 








1 


8 = 0, e < 


N-i 


1 





kl 


8 < 0, e > 


1 








1 


8 < 0, e = 


1 








1 


8 < 0, e < 


N-i 


1 





£ 



TABLE II: Number n v (8,e) of up interfacial spins of class 
77 = 2, 3, 4 for the column configuration types defined in Table 
Q The last column shows the total number of up interfacial 
spins associated to each type. 



This completely characterizes the population of each in- 
terfacial class in the involved column. Consider, for in- 
stance, a case 5, e > 0, as in the column between x — 4 
and x = 5 in Fig. ^ This column contains S + e in- 
terfacial spins, of which 8 + e — 2 are of class 2 and the 
other two spins are of class 3. On the other hand, this 
configuration involves 5 + e + 1 broken bonds, of which 
5 + 1 belong to up spins in this column; following our 
convention above, the other e broken bonds may be asso- 
ciated to the up interfacial spins in the column between 
x + 1 and x + 2. In order to go further in the analytical 
solution of the problem, let us neglect column-column 
correlations by assuming that broken bonds in the inter- 
facial column are to be associated to interfacial spins in 
this column. For S, e > 0, one has S + e — 2 broken bonds 
associated to interfacial spins of class 2, and three bro- 
ken bonds associated to spins of class 3. Let us assume 
also that, in general, the probability of a given interfa- 
cial column configuration is proportional to the product 
of probabilities of each of the broken bonds which form 
it. We have that 



P(S > 0,e > 0) : 



i{A^X|+- 2 X|}, 



= e 7W) .E 



where Q is a normalization factor, and A(</>) 
Equivalently, for 6 > and e < one may write that 



P{6 > 0,e < 0) 



1 



{A- 5 



A 3 



} 



where a = min(|<5|, |e|) and A = max(|<5|, |e|). Table [I] 
shows all of the possible interfacial column configura- 
tions, together with their statistical weights. Notice that, 
as one may conclude from Tableland eq. the proba- 
bilities P(6, e) converge as p — > to the SOS equilibrium 
value XP(S)P(e). 

The densities H v (T,p) may be written as an average 



5 



over all possible interfacial column configurations: 

oo 

U V (T )P )= £ 7T v (6,e)P(6,e). (9) 

5,e— — oo 

Here, 7r^((5, e) is the probability of finding an up spin of 
class 77 € [2, 4] in an interfacial column characterized by 
the pair (<5, e). In general, tt v (S, e) — n ri (d,e)/N(d,e), 
where n v (6, e) is the number of up spins of class rj in 
an interfacial column characterized by (5, e), and N(5, e) 
is the total number of up interfacial spins associated to 
this column. Table ITT1 shows rii(5,e) and N(5,e) for all 
possible configurations. 

The densities (0 will depend on the average inter- 
face slope, tan </>, through the Lagrangian multiplier 7(1/)). 
This dependence is inherited in general by the average 
weight of a broken bond in the noncquilibrium regime, 
X(T,p), see eq. (|SJ), and it makes the explicit calcula- 
tion of the noncquilibrium surface tension unfeasible, sec 
§ III. a. Therefore, further simplifications are needed. We 
shall assume that the densities II,, (T, p) in eq. corre- 
spond to the case tan(/> = 0, i.e., A(</>) = e 7< ^) = 1. In 
fact, the underlying lattice anisotropy implies that the in- 
terface tends in general to orientate parallel to any of the 
principal axis, which has a lower energy cost. Therefore, 
for regions of the parameter space (T, p) where such ten- 
dency is strong, it is justified to particularize the popula- 
tions H v (T,p) to the case tan0 = 0. On the other hand, 
the parameter-space regions in which that tendency is 
weak are characterized by an effective isotropy, so that 
particularizing to a given orientation, e.g., tan <f> = is 
valid. This approximation amounts to assume that the 
relevant orientation dependence entering the definition of 
the nonequilibrium surface tension, eq. Q), comes from 
the dependence on A(<^) that appears in the probabil- 
ity P{8) of a step of size 5 in the interface, see eq. 
A higher order, iterative procedure to take into account 
the ^-dependence of X(T,p > 0) would consist in: (i) 
calculate X^ = X(T,p) using the tan(/) = simplifi- 
cation, (ii) use to compute A(°)((/)) as a function 
of <j) from eq. JHJ), (Hi) replace A^(cf>) in the general, 
A(</>)-dependent expression for X(T, p), and use the so- 
defined X^- 1 ' in Q to compute an improved approxima- 
tion to the nonequilibrium surface tension. We do not 
expect this complex procedure to produce a significative 
improvement of our approximation (while the resulting 
formulas are much more involved). 

Assuming A(0) = 1, and using 7r^(<5, e) as given in Ta- 
bic |nl we obtain from eq.® that 



The same method leads to 



n. 



2X a 



X : 



+ 



■ln(l-X 2 



X* 



(1 - x 2 y 1 - x 2 
x 2 x 3 



2 + ^ln(l-X 2 ) 



X 



1-Xo 



X 2 — Xr, 



1-Xf 



(10) 



Xn 



xl 



X 2 — Xi 



In 



l-X 2 
1-Xf 



1-X$ 



X2 
xl 



n. 



X^Xo 



Q (X 2 - XI 



\ 



111 



1-Xf 



X2 (1 - x 2 y 



In ( 1 - Xi 



(11) 

and the normalization condition gives Ii3(X',p) = 1 — 
H 2 (T,p) — Ii4(T, p). The factor Q follows from the nor- 
malization of P(S, e) for A((/>) = 1 as 



s = 



2Xj 
L — X 2 

(1- 



Xo 



1-X 2 



Xn)Xi 



(x 2 -xi)(i-xi) 



x^x 2 
x 2 - xl 



Xj,X 2 

~ ~xi 



X, 



In equilibrium, p = 0, where X v (T,p) — > X(T,p = 0) = 
e _2/3 , the above expressions for H v (T,p) reduce to the 
known SOS equilibrium results. [Isj 



IV. SOME RESULTS 

We summarize in this section some main results that 
follow from the above for the properties of the nonequilib- 
rium interface, and compare our predictions with Monte 
Carlo simulation data. 



Microscopic Structure 

The basic SOS hypothesis is that the probability of a 
step of size S in the nonequilibrium interface is given by 
©. That is, P(S) is an exponentially-decaying function 
of \6\ controlled by a typical scale which, for tan0 = 0, 
is 6 — \hiX(T, p)!" 1 with X given in (JSJ . In order to 
check this assumption on the microscopic structure of 
the interface, we performed Monte Carlo simulations of 
the system in § II. This evolved with time in the com- 
puter starting with two stripes of the same width of up 
and down spins, respectively, separated by a flat inter- 
face. We suppressed bulk dynamics in these simulations, 
i.e., bulk, class-1 spins remained frozen to prevent the 
nucleation of droplets in the bulk to interfere the inter- 
face dynamics. This turns out to simplify notably the 
analysis while it does not modify essentially the interface 
structure except close to the critical temperature, where 
fluctuation of all sizes occur. In fact, we carefully checked 
the validity of this assertion by computer simulations. 

The interface was thus observed to eventually reach 
a steady state, in which we measured the inter- 
face microscopic structure as a time average of S = 
{Si, i = 1, . . . , L x } — taking into account some (small) 
correlations observed between neighboring steps. Fig. |3 
depicts P(S) as obtained for a (relatively large) system 
at low temperature for an average interface slope such 
that tan(/> = 0. The figure shows also our theoretical 
prediction, eq. ©, revealing an excellent agreement for 



6 




FIG. 2: The symbols are Monte Carlo results for the prob- 
ability of step 8 for a system of size L x x L y — 256 x 128 
at temperature T = 0.3Tb, with tan<j!> = and different val- 
ues of the nonequilibrium perturbation, namely, p — 0, 0.01, 
0.02, 0.03, 0.04 and 0.05. from bottom to top. Solid lines are 
the corresponding teoretical prediction JHJ. For the shake of 
clarity, the curves are shifted by a factor 2 % , i £ [0,5] in the 
vertical direction, where i = 100 x p. 



all values of p. It is noticeable that the typical step scale 
8 increases as p increases at fixed T, i.e., the nonequilib- 
rium noise tends to amplify the interface fluctuations, as 
one should have probably expected. 

Our data is good enough to provide an estimate for the 
second central moment of the step distribution, which 
measures the interface width, 
Our theoretical prediction is 



! (7» = (S 2 ) - tan 2 



w 2 (T,p) 



XA 
2XA 



1 



1 



(l-XA) 2 (A-X) 2 
2X 



(l-XA) 3 (A-X) 3 



— tan^ 



where A((f>), z((f>) and X(T,p) are given by eqs. J5J, © 
and ©, respectively. Fig. [3] compares this with Monte 
Carlo estimates. The most noticeable fact here is that 
w 2 (T, p) extrapolates towards a non-zero value in the low 
temperature limit for any p > 0, contrary to the case 
of an equilibrium interface, which is completely flat at 
zero temperature. That is, nonequilibrium fluctuations 
imply a rough interface even at zero temperature. The 
low temperature roughness may be estimated by realizing 
that, for T — > and moderate values of p, Il2(T, p) ~ 1 
from eq. tflOpt . i.e., almost all interfacial up spins belong 
to class 2. Therefore, X(T,p) -> X 2 (T = 0,p) = </p in 
this limit, so that one has for tanoi = that 



w 2 (T = 0,p) 



(13) 



The inset in Fig. depicts this behavior which is in full 
agreement with our Mote Carlo values. The only signifi- 
cant differences we found between theory and simulations 




TAi\ 



FIG. 3: The symbols are Monte Carlo results for the interfa- 
cial width w 2 (T,p) as a function of temperature for a system 
of size L x x L v — 256 x 128, tan</> = 0, and, from bottom to 
top, p — 0, 0.005, 0.01 and 0.02. Errorbars are smaller that 
the symbol sizes. Solid lines are the corresponding theoretical 
prediction. Notice the non-zero interfacial width in the low 
temperature limit for the nonequilibrium system (p > 0) . The 
inset shows the SOS zero-temperature limit for the interfacial 
width, w 2 (T = 0,p), see eq. 1131 , as a function of p. 



are at high enough temperatures, as illustrated in Fig. [3] 
This is to be attributed to step-step correlations as the 
critical temperature is approached. 



Macroscopic Behavior 



A principal interface macroscopic property is the sur- 
face tension. The prediction Q is illustrated in Fig. 0]for 
6 (12) 6 = 0. It is remarkable the essential difference occurring 
at low T between the equilibrium and nonequilibrium 
cases. 

The surface tension cr ne (T, p > 0; 4> = 0) exhibits non- 
monotonous behavior as a function of T, with a maximum 
at a temperature which depends on the intensity of the 
nonequilibrium perturbation, T max (p). For T < T ma x(p), 
CT nc decreases as one cools the system. This anoma- 
lous behavior turns out to play a fundamental role in 
understanding the exit from metastable states in this 
system. [HI ll'J and it is likely it may help the under- 
standing of the low temperature behavior in other com- 
plex systems concerning nucleation and growth processes. 



We devised the following indirect method to check 
our predictions for the nonequilibrium surface tension. 
The system defined in § II happens to exhibit long-lived 
metastable states in the ordered phase when subject to 
a small negative external magnetic field. [I^ITsT ] The exit 
from this state is a highly inhomogeneous process that 
proceeds via the nucleation and growth of one or sev- 
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T/T„ 



FIG. 4: Main graph: Theoretical prediction for the surface 
tension as a function of temperature for, from top to bottom, 



0, 1CT 



10" 



10" 



10~ 3 , 5 x 10" 3 , 



3 x 10~ and 4 x 10~ . Notice that, for any — even small — 
p > 0, the surface tension behaves non-monotonously, con- 
trary to the equilibrium case. Inset: The effective interface 
temperature, T ef ■ , as defined in the main text, as a function 
of T for the same values of p than in the main graph. No- 
tice that T^{T,p > 0) strongly deviates from T in the low 
temperature regime. 



eral droplets of the stable phase within the metastable 
sea. Droplet nucleation is controlled by the compe- 
tition between the surface tension, which hinders the 
droplet growth, and the bulk "free energy" , which favours 
it. Consequently, small droplets — having a large sur- 
face/volume ratio — tend to shrink, while the larger ones 
tend to grow. The critical droplet size, R c (T,p), sep- 
arates these two regimes. Following further the trend 
in equilibrium theory. [l^| one may assume that an effec- 
tive macroscopic potential controls the escape from the 
metastable state for < p < 1, and that 



Rc(T,p) 



(d-l)a(T,p) 
2m s (T,p)\h\ ' 



(14) 



Here, d is the system dimensionality, h is the applied 
magnetic field, a(T,p) stands for the zero-field surface 
tension along one of the lattice axis, and m s (T, p) is the 
spontaneous (positive) magnetization for h — 0. The lat- 
ter may be approximated by mean field-theory, 0] and 
the surface tension may then be obtained from a Monte 
Carlo estimate of R c (T,p). 

In order to perform this computation, consider our sys- 
tem in a square lattice with periodic boundary conditions 
along both the x and y directions. The Hamiltonian is 
now Ti' = H — hj^i s ii with Ti given in Q and h < 0. 
All the spins are initially up, except for a square droplet 
of down spins of side 2R which represents the stable 
phase. This is let to evolve according to (0- The state 
is highly unstable, so that any subcritical initial cluster, 
i.e., R < i? c (T, p), will very quickly shrink, while a super- 
critical one, R > R c (T,p), will rapidly grow to cover the 



FIG. 5: Critical droplet size, 1Z C , as a function of temperature 
for a system of size L = 53, with periodic boundary condi- 
tions, subject to a magnetic field h = —0.1, and, from top to 
bottom, p = 0, 0.001, 0.005 and 0.01. The symbols are Monte 
Carlo results obtained as an average over N exp = 1000 inde- 
pendent "experiments". The solid curves correspond to the 
theoretical prediction l|140 . For the shake of clarity, results for 
the n-th value of p, n — 1, . . . , 4 (using the above indicated 
order) have been shifted by (1 — n) units along the y axis. 
The inset shows Monte Carlo results for the probability that 
a droplet of radius R is supercritical, P spc (R), as a function 
of R for a system of size L = 53 at T — 0.4Tb, p = 0, and 
h — —0.1. This corresponds to 10 3 independent experiments 
for each value of R. In all cases, error bars are smaller than 
the symbol sizes. 



whole system. Since our dynamics is stochastic, we define 
the probability that a droplet of size R is supercritical, 
P spc (R). This is measured in practice by simply repeating 
many times the simulation and counting the number of 
times that the initial droplet grows to cover the system. 
The critical droplet size is defined by P spc (Rc) = 0.5. As 
observed in the inset of Fig. [5] P spc (R) sharply goes from 
to 1, which allows a relatively accurate estimate of R c . 

Fig. [5] compares our Monte Carlo results for R c (T,p) 
with the theoretical prediction from (|14|l using a(T,p) — 
a ne (T,p; <p — 0) as given by eq. (0. The agreement is 
rather good, and we confirm that R c (T,p) behaves non- 
monotonously with T in the nonequilibrium regime. 

The temperature dependence of er nc (T, p; cf> = 0) may 
be used to compute the phase diagram of the model, 
Tc(p). 24] In equilibrium, the interface free energy ap- 
proaches zero as T — > To; for T > To there is no surface 
tension because there exist only one disordered phase. [2^ 
Therefore, if as assumed in this paper a nc captures the 
macroscopic properties of the nonequilibrium interface, 
we may identify Tc{p) as the temperature (other than 
T = 0) for which cr no (T,p,0 = 0) = 0. This is done in 
Fig. El In particular, we may ask about the nonequilib- 
rium parameter p c above which no ordered phase exists 
at low T. For p > and T — > 0, a ne (T,p; = 0)- a(p)T 
(see Fig. 0] and eq. Q). The slope a(p) decreases 
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o.s 



Ordered phase 



Disordered phase 




P C =(V2-1) : 



we conclude that 



0.15 



0.2 



FIG. 6: Phase diagram of the model as obtained from the 
nonequilibrium surface tension. Notice the reentrant behavior 
of T c {p) for p c < p < p*, with p c = {y/2 - l) 2 and p* w 
0.18625. Inset: a nc (<j) = 0) as a function of T for p = 0.175 > 
p c - Notice the negative slope at low-T and the intermediate 
temperature regime where <r nc is positive. 



monotonously with p, and the condition a(p c ) = sig- 
nals the onset of disorder at low temperature. This yields 
p c = (V2— l) 2 s=s 0.1716, in excellent agreement with pre- 
vious Monte Carlo simulations. |Toj For p > p c one does 
not expect low-T order. However, the non-monotonous 
temperature dependence of a nc (see inset in Fig. HJ) 
involves the emergence of an intermediate-T region for 
p c < p < p* « 0.18625 where order sets in, as opposed to 
the low-T and high-T disordered phases, see Fig. EJ This 
reentrant behavior of Tq (p) is similar in spirit to the one 
reported in Ref. 13] for the spinodal held characterizing 
the limit of metastability in this model, and it is reminis- 
cent of the reentrant phase dia gra m observed in systems 
subject to multiplicative noise. [25j 

In order to gain some intuition on the physical origin of 
the anomalous low-temperature behavior of the nonequi- 
librium surface tension, let us write the statistical weight 
of an interfacial broken bond as 



X(T,p) = exp 



(15) 



which defines an interface effective temperature. The 



function TjjjP (T,p) = l/P^' is illustrated in the inset of 
Fig. H One hrst realizes that T^ } (T,p > 0) > T for 
any T £ [0, To]. That is, the nonequilibrium interface 
endures an effective temperature larger than the thermo- 
dynamic one. On the other hand, one identifies two dif- 
ferent regimes in the inset of Fig. 0]at given p. At high-T, 
where thermal fluctuations dominate over the nonequi- 
librium noise, is proportional to T. However, at low 
enough T the nonequilibrium noise dominates; in this 
case T e £ deviates from T and tends to a saturating, con- 
stant value which depends on p. Following the method 
above to obtain the zero-temperature interfacial width, 



Km TV(T,p)n-— >0. 
t^o 01 lnp 



(16) 



It is also remarkable that the onset of the deviation of 
T j from T for a given p coincides with the maximum 
observed for the nonequilibrium surface tension; see Fig. 
H In fact, since T c ^ is a small constant for < T 
Tc(p) and p <C p c , we may expand a ne at low-T to obtain 
cr„ c ~ (T/TjfVc + 0(TX). On the other hand, for < 
T < T c (p) and the same p < p c , the quotient T/T^ = 
1 — a is a constant, with < a « 1. Using a as small 
parameter now, a high-T expansion of er nc yields er nc ~ 
(T/T e f )<7 e + O(a). In both limits the corrections to the 
asymptotic behavior er nc ~ (T/T^j')a e are small. Now, 
since T c ( f :) ~ constant and cr c ~ 2(1 - Te~ 2 / T ) as T -> 
0, one expects a ne to be an increasing function of T in 
this limit. On the other hand, T/T^ 1 ~ constant in 
the high-T limit, so er no depends on T as a e does, i.e. 
Cne decreases with T for high enough T. Therefore one 
would expect a non-monotonous T-dependence of a ne , 
with a maximum at T max (jp), which roughly coincides 
with the crossover observed in T^(T,p). In this way, 
the anomalous T-dependence of a nc can be traced back 
to the crossover between a T-dominatcd, high-T regime 



and a p-dominated, low-T region, as captured by T 



(i) 

cf • 



Droplet Shape 



The droplet shape is controlled by the need to min- 
imize the total surface tension at constant droplet vol- 
ume. For isotropic systems, this implies spherical shape. 
In our case, however, the surface tension depends on the 
orientation of the interface with respect to a privileged 
axis, a(4>). Consequently, the shape adjusts itself to take 
advantage of the low free energy cost of certain inter- 
face orientations, which produces droplets with a crystal- 
like appearance which depends on temperature and other 
parameters. |26J We apply next the Wulff construction 
|27| to obtain information concerning the nonequilibrium 
droplet shape. 

The method essentially consists in considering the po- 
lar curve a(4>) 1 <f> € [0, 2ir], and drawing through its points 
a line perpendicular to the radius. The interior envelope 
to these lines determines the droplet radial function in 
polar coordinates, R(9). More specifically, one may write 
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FIG. 7: (a) Shape of a droplet, as obtained from the Wulff 
construction, for p — (equilibrium) at, following to the cen- 
tre, T/T = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 and 0.9. For 
the shake of clarity we have rescaled the droplet according to 
its temperature, (b) The same as in (a), but for the nonequi- 
librium model with p = 0.01. 



parametrically: [5^ 

me) = i? [z 2 W>) + ?/W /2 , 
x{4>) 



(17) 



<y\cp) cos <p — sin < 



2/(0) 



tan 9 



cr{<p) sin© — — cos< 

dcp 

y{4>) 
x&y 



where R is a fixed length scale, and ct(0) is the sur- 
face tension. We approximate the latter by Q. This 
is needed only for the angular interval <f> S [0, 7r/4], 
since one may extend then to the whole circumference 
by straightforward symmetry considerations. The result 
is singular for angles 4> — (2n+l)7r/4, n = 0, . . . , 3, which 
gives rise to angular intervals around 9 = (2n + l)7r/4, 
n = 0, . . . , 3, where R{9) is not defined. Therefore, one 
considers the analytical continuation r(9) such that, in 
particular, dr/d9 — at 9 = n/A as required by sym- 
metry. This, together with continuity and analyticity, 
leads to a second order polynomial and its coefficients, 
r{9) = a9 2 + b9 + c. 

Fig. [7| illustrates the result. In equilibrium, p = 0, 
the droplet tends to become squared as T — ► due to 
the underlying lattice anisotropy, while it recovers the 
(isotropic) spherical shape for T > 0.5Xb- In nonequilib- 
rium, the droplet adopts a shape which is intermediate 
between a circle and a square. This is again due to the 
fact that the temperature T?' that the interface feels 
does not go to zero as T — ► for any p > 0. 

A more quantitative description is provided by the 
droplet form factor, fld(T,p). This is defined via the 
equality V = ild(T,p)lZ d , where 1Z = R(9 = 0) is a mea- 
sure of the droplet radius, and V is the droplet volume. 
For a two-dimensional system, 




(18) 



T/TV 



FIG. 8: The form factor Q(T, p) as a function of tempera- 
ture for, from top to bottom, p = 0, 0.001, 0.005, 0.01, 0.02 
and 0.03. The top (bottom) line corresponds to the squared 
(circular) droplet. 



The square and circular droplets are characterized by 
£1 = 4 and tt, respectively. Fig. [S] shows fl(T,p) as a 
function of T for different values of p. This clearly demon- 
strates that (T, p) goes to 4 (squared shape) in the low- 
T limit for p = 0, but the tendency is towards a smaller 
value for any p > 0. That is, unlike in equilibrium at low 
enough temperature, no facets are expected in a nonequi- 
librium droplet. 



CONCLUSIONS 

This paper deals with the influence of a simple 
nonequilibrium condition, which may ideally represent 
the situation in a class of disordered systems, [lfj on 
the microscopic and macroscopic properties of a com- 
plex interface. It is assumed that the probability of a 
discrete change by d at the (nonequilibrium) interface is 
proportional to X' 5 !, where X = X{T,p) is the statisti- 
cal weight of a interfacial broken bond for temperature 
T and nonequilibrium perturbation p. This is expressed 
X = J^r) n»;^?? m terms of the probability II,, that an 
interfacial broken bond ends at a spin surrounded by cer- 
tain degree, 77, of local order. 

Our main hypothesis consists in assuming that the 
nonequilibrium system is attempting to minimize a sur- 
face and, consequently, one may translate here the equi- 
librium formalism. It is also assumed that, at least for 
the model studied in this paper, the specific nonequilib- 
rium mechanism simply adds to the basic thermal mech- 
anism in such a way that it may be incorporated in a 
non-perturbative manner to the microscopic parameter 
X(T,p). Therefore, this contains all the information con- 
cerning the effect of the nonequilibrium perturbation p 
on the interface, and it follows that a SOS theory based 
on X(T,p) yields the micro- and macroscopic behavior of 
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the nonequilibrium interface. In this way, one may finally 
obtain an explicit expression for the relevant nonequilib- 
rium surface tension, a ne (T, p) . 

Regarding the microscopic interface structure, the 
nonequilibrium noise turns out to enhance interfacial 
fluctuations. In particular, the typical scale for interfa- 
cial fluctuations increases with p. It is also demonstrated 
that the nonequilibrium interface remains rough in the 
zero-T limit, contrary to the equilibrium case. These 
are theoretical predictions in full agreement with Monte 
Carlo simulations. 

Regarding macroscopic behavior, o~ nc (T,p) exhibits 
anomalous behavior at low T (for any p > 0). In par- 
ticular, cr no is a non- monotonous function of T with a 
maximum at T — T max (p), and er no decreases as the 
system is cooled further below T max (p). This counter- 
intuitive prediction is also confirmed indirectly by Monte 
Carlo simulations. That is, we estimated numerically the 
critical droplet size, R c (T, p) , which is expected to be pro- 
portional to the surface tension. |l2| as the system exits 
from a metastable state. Some intuition on the origin of 
this anomaly is obtained by defining an interface effective 
temperature which importantly deviates from T. In this 
way, the non-monotonous T-dependence of o~ ne can be 
related to a crossover between two different temperature 
regimes: a low-T region dominated by the nonequilib- 
rium noise, where o~ nc cx T, and a high-T regime dom- 
inated by thermal fluctuations, where o~ nc cx o~ c . The 
shape of the nonequilibrium droplet as obtained by a 
Wulff construction also reflects the anomaly of a nc - We 
find, in particular, that droplets at very low temperature 
tend to minimize more their surface under the nonequi- 
librium condition. 

These details are essential to nucleation theory. There- 
fore, we expect that the anomalous low-temperature be- 
havior of the nonequilibrium surface tension described 
above may be relevant to many physical processes such 
as the ones mentioned in § I. The possible utility of our 
results here will be addressed in a forthcoming paper con- 
cerning the relaxation of a nonequilibrium system from 
a metastable state, [l^ 

Finally, the results in this paper are explicitely ob- 
tained for a square lattice. Some caution should be used 
before generalizing, since there are examples when the 
shape and properties of a nonequilibrium interface de- 
pend strongly on the geometry of the host lattice. |2flj 
However, we believe that the phenomenology here de- 
scribed should hold for more lattice geometries other than 
square, provided that the (nonequilibrium) interface en- 
dures an effective temperature with the same qualitative 
properties than the one discussed above, i.e. T^jP satu- 
rates to a constant, p-dependent value as T — > and is 
proportional to T at high enough temperature. 
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